High Curie temperatures in (Ga,Mn)N from Mn clustering 
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Abstract 

The effect of microscopic Mn cluster distribution on the Curie temperature (Tc) of (Ga,Mn)N 
is studied using density-functional calculations together with the mean field approximation. We 
find that the calculated Tc depends crucially on the microscopic cluster distribution, which can 
explain the abnormally large variations in experimental Tc values from a few K to well above room 
temperature. The partially dimerized Mn2-Mni distribution is found to give the highest Tc > 500 
K, and in general, the presence of the Mn2 dimer has a tendency to enhance Tc- The lowest Tc 
values close to zero are obtained for the Mn4-Mniand Mn4-Mn3 distributions. 
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The discovery of ferromagnetism in (In,Mn)As [lj sprouted a family of III-V based diluted 
magnetic semiconductors. According to the mean-field Zener theory (Ga,Mn)N should have 
the highest^ Curie temperature Tc (— 400 K) among these prospective III-V spintronics 
]. However, the experimental Tc values for (Ga,Mn)N range from 10 to 940 
. First principles calculations combined with the mean field approximation have 
given rather high Tq values of about 270-350 K 

Baa 

, while more sophisticated statistical 
approaches includin g pe rcolation effects and/or magnetic fluctuations give significantly lower 
Curie temperatures 10, The substantial variation of the experimental T c values 




is suggested to be due to dustering of Mu atoms Bill and the formation of giant magnetic 
moments at the Mn clusters . Consistently with this, it is found tha t it is energetically 
favourable for substitutional Mn atoms to cluster around N atoms 15. 101. Also, Mn clusters 
are observed directly at high Mn concentrations 
ca 
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la 19. Although in many theoretical 



21 



10 



20 



culations Mn clustering is generally found to decrease T^'s in III-V materials ?1 
l^ . we show in this Letter by using first principles calculations and the mean field 
theory that Mn clustering may also act beneficially by increasing Tc- More generally, we 
find that the stability of the ferromagnetic order depends crucially on the microscopic cluster 
distribution: small Mn 2 clusters may increase Tc whereas Mn 4 clusters may reduce Tc down 
to a few Kelvins (by an Mnj cluster with i = 2,3,4 we mean i Mn atoms surrounding the 
same central neigbouring N atom). 

The presence of Mn clusters in (Ga,Mn)N is unavoidable at the Mn concentrations that 
give high Tc values. This is obvious from Table H] where Mn cluster distributions are given 
when the Ga atoms are substituted randomly by the Mn atoms at various relevant Mn 
concentrations. Already at x — 0.02 (counting single Mn atoms as Mni clusters) the Mn 2 
cluster portion exceeds 10 % and grows fast up to 40 % with increasing x. Also the Mn 3 
portion increases markedly with increasing x whereas the Mn 4 cluster portion remains less 
than 0.5 % even for highest Mn concentration of x = 0.14. To investigate the microscopic 
cluster distribution effects on the Curie temperature we perform spin-polarized total energy 
density-functional calculations using the projector augmented wave method as implemented 
in the VASP code j^. Our supercells contain two Mnj clusters (i = 1,2,3,4) at maximum 
separation, thus representing a uniform distribution of clusters. The exchange-correlation is 
approximated as follows. In the first approach we use the generalized gradient approximation 
(GGA-PW91). In the second approach we use the local spin density approximation plus 
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XI (L + XI) with U chosen as 3 eV and used only for Mn. The second approach tries to 
compensate the self-interaction which is present in the local density approximation and 
which is known to be substantial in (Ga,Mn)N j^j. The plane-waves cut-off is 425 eV. We 
use supercells of 48, 72, or 108 atoms (3x2x2,3x3x2, or 3x3x3 primitive wurzite unit 
cells, respectively) and the 4x6x4, 4x4x4 or 4x4x3 fc-meshes for the Brillouin zone 
sampling, respectively. We use throughout the same lattice parameters a = 3.217 A and 
c : a = 1.631 obtained from geometry optimization of the wurzite GaN crystal. The atomic 
positions are fully relaxed in the 48 and 72 atom supercells, and in the 108 atom supercells 



the atomic positions are fixed to the respective ones of the smal 
Tc values are roughly estimated from the mean-field expression 



er relaxed supercells. The 
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where iV is the number of Mn clusters in the supercell fixed to 2, and A is the difference in 
supercell total energy between the anti- par allel and parallel alignments. 

First we find, in agreement with Ref. 16j, that the Mnj cluster (i = 2,3,4) is energetically 
the most stable configuration of i Mn atoms. This suggests that the Mn 2 , Mn 3 and Mn 4 
cluster portions may in reality be even larger than the ones based on random substitution 
given in Table HI We also find in agreement with Refs Q) that large stable magnetic 
moments are formed at the clusters, and therefore intercluster interactions determine the 
ferromagnetic order. 

We now address the clustering effects on the Curie temperature within the GGA and the 
mean field approximation. In the case of the (nearly) uniform distribution of the single Mn 
atoms, when the Mn concentration x increases from 0.056 to 0.083 (and the shortest Mn 
- Mn distance decreases from 7.65 to 6.15 A) the calculated Tc increases sharply from 39 



to 355 K, as seen in the first two rows of Table ITU qualitatively in agreement with Ref. jjj]. 
The increased magnetic Mn - Mn coupling is reflected in the density of states (DOS) (given 
in Fig. ^ (a)) where the width of the midgap Mn d type majority spin DOS increases from 
0.5 to 1 eV (see the inset in Fig. ^ (a)). Next we allow partial dimerization at the same 
concentration of x = 0.083. Two thirds of the Mn atoms are allowed to form Mn 2 clusters 
in the plane perpendicular to the c axis. This causes a further increase of Tc up to 514 
K given in Table |H] In Fig. Q (b) for x = 0.083 the partial dimerization is seen to cause 
further broadening of the midgap majority spin DOS into a 1.7 eV wide joined structure 
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with an increased p - d hybridization. This atomic configuration offers the highest Tc we 
have found. In more detail, the DOS of a uniform Mn 2 distribution for x = 0.056 is given in 
the inset in Fig. 0(b), and comparison with this inset shows that the structure in the DOS 
is mainly determined from the Mn2 DOS. 

To see whether clustering from Mn2 to Mn 3 could raise the Curie temperature further, 
above 514 K, we add a third atom to the Mn 2 cluster in the plane perpendicular to the c 
axis. However, in this case the clustering decreases Tc to 294 K in agreement with earlier 
[2, 10, 20, 21, 22! . The midgap DOS splits into four distinct subbands (with three 
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gaps) still with a clear p - d hybridization. When the distance between the Mn3 cluster and 
the single Mn atom is increased (the distance between the central N atom of the cluster and 
the single Mn atom grows from 5.90 to 9.08 A) Tc decreases further to 128 K (see Table ITT|) . 
as expected. 

Finally, we complete the Mn 3 cluster to a symmetric Mn 4 tetramer cluster and find 
that Tc diminishes down to 27 K (Table ITTj) . The corresponding DOS shown in Fig. (a) 
consists of the narrow peak of about 0.5 eV at the Fermi level and of several other separated 
peaks. The total width of the original majority spin midgap DOS of 0.5 eV has grown to 
2.1 eV which exceeds the calculated gap. The upper part overlaps with the minority spin 
conduction band edge. The reason for the low Tc may be understood by comparing with 
the DOS of a uniform Mn 4 distribution (x = 0.093) shown in Fig. (b). It is immediately 
obvious that the essential p - d hybridization is completely missing from the narrow peak at 
the Fermi level. The central narrow subband is of purely Mn d type because the N p orbitals 
are mainly used for the Mn 4 intracluster bonding which makes the magnetic coupling of the 
Mn 4 clusters to other Mn clusters (in this case to Mni) weak. Summarizing, we find that 
the presence of the Mn 4 clusters leads to a low Curie temperature. This explains also why 
in the Mn 4 - Mn 3 case Tc drops to 4 K although the Mn concentration is as high as x = 0.13 
(Table . On the other hand, the presence of the Mn 2 clusters has a tendency to increase 
T c . 

The L + U results differ from the GGA ones mainly in shifting the midgap Mn d type 
DOS towards the valence band maximum by 0.5 eV. In the partial dimerization the Mn d 
type DOS merges with the top of the valence band causing a growth of Tc up to 618 K 
(Table ITT]) . This behaviour begins to resemble that in (Ga,Mn)As and might be assigned 
to the p - d type hole mediated ferromagnetic coupling. In the Mn3-Mni, and especially 
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in the Mn 4 -Mn 1 case, the p - d type empty peak splits off decreasing T c to 453 and 39 K 



221. All in all it seems that the 



respectively ( Table ITT]), similarly as for (Ga,Mn)As in Ref. 
previous conclusions on clustering effects on the Curie temperature hold also in the case 
when further corrections to the electronic correlation beyond L + U are included. 



The mean fie 



d approximation used in this work is known to overestimate Curie tem- 



peratures i.e. the dependence of Tc on the total energy difference A may 

not be linear as in Eq. (0). Nevertheless, high A will imply high Tc, and thus the general 
trends presented for different microscopic cluster configurations should hold. It is interest- 
ing to remark that a recent study considering dimerization effects beyond the mean field 
approximation leads to a similar conclusion as the present work. 

In conclusion, we find that the calculated Curie temperature depends crucially on the 
microscopic cluster distribution. The partially dimerized Mn 2 -Mn! distribution is found to 
give the highest Tc- This suggests that the presence of the Mn 2 dimer provides a mechanism 
to enhance Tc- On the other hand, low Tc values are obtained for the Mn4-Mni and M114- 
Mn3 distributions. Thus we have shown that different cluster distributions may explain the 
abnormally large variations in experimental Curie temperatures. 
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TABLE U] Spin-flip energies (A) and Curie temperatures (Tq) for two clusters in the 
supercell. The fourth and fifth columns are GGA calculations, the sixth and seventh columns 
h+U calculations (denoted by the index "+[/"). d is the minimum distance between the 
Mn atoms belonging to the different centers. 

TABLE [H] Spin-flip energies (A) and Curie temperatures (Tc) for two clusters in the 
supercell. The fourth and fifth columns are GGA calculations, the sixth and seventh columns 
h+U calculations (denoted by the index "+[/"). d is the minimum distance between the 
Mn atoms belonging to the different centers. 

FIG[0 Densities of states together with Mn d and N p orbital projections for most likely 
cluster arrangements in (Ga,Mn)N. (a) Uniform Mn distribution for x = 0.056, in the inset 
for x = 0.083. (b) Partially dimerized Mn2 - Mni for x = 0.083, in the inset a uniform Mn2 
distribution for x = 0.056. 

FIG Densities of states for (Ga,Mn)N. (a) Mn 4 - Mm for x = 0.111. (b) Uniform Mn 4 
distribution for x = 0.093. For further information see the caption of Fig. ^ 
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TABLE I: Mn cluster distribution as a function of x when the Ga atoms are substituted randomly 
by the Mn atoms. The portions are given in units of % (counting single Mn atoms as Mni clusters). 

x Mni Mn 2 Mn 3 Mn 4 

0.02 89 11 

0.06 71 27 2 

0.10 58 37 5 

0.14 49 42 9 
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TABLE II: Spin-flip energies (A) and Curie temperatures (Tq) for two clusters in the supercell. The 
fourth and fifth columns are GGA calculations, the sixth and seventh columns L+J7 calculations 
(denoted by the index "+[/"). d is the minimum distance between the Mn atoms belonging to the 

different centers. 

System x (%) d (A) A (meV) T c (K) A +u T+ u 



Mni - 


Mni 




5.6 


7.65 


10 


39 


1 
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Mni - 


Mni 
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8.3 


6.15 


92 


355 






Mn 2 - 


Mni 
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8.3 


6.13 


133 


514 


160 


618 


Mn 3 - 


Mni 
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11.1 


6.13 


76 


294 


117 


453 


Mn 3 - 


Mni 
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7.4 


9.28 


33 


128 






Mn 4 - 


Mni 


c 


9.3 


7.66 
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27 


10 


39 


Mn 4 - 


Mn 3 
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13.0 


6.19 


1 
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a 72 atoms/supercell (SC), b 48 atoms/ SC, c 108 atoms/SC 
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